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Abstract 

The vacuum polarization tensor in strong external magnetic fields has been 
evaluated numerically for various strengths of magnetic fields and momenta 
of photons under the threshold of the pair creation. The fitting formula 
has been obtained which reproduces the calculated results within 10 % of 
error. The proper time method is employed further to obtain the retarded 
polarization tensor for finite temperature plasmas. 
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I. INTRODUCTION 



The strong magnetic field is attracting attentions of astropliysicists tliese days. It fias 
been known tliat tlie magnetized vacuum sliows interesting features as tfie magnetic field 
strength exceeds a critical value = ml/e ~ 4 x lO^^G [|l],H. Since this value is so large 
even in the universe compared, for example, with the canonical magnetic field of lO^^G for 
a pulsar, that it was supposed that this was a subject of academic interest only. This has 
been changing drastically recently. 

Some observations suggest the existence of neutron stars with a magnetic field far 
greater (~ lO^^G) than the canonical one for the observed pulsars, and they are called a 
magnetar as dubbed by Duncun and Thompson 0]. As the reality of very large magnetic 
fields looms, some researchers speculated further that some other peculiar extraterrestrial 
phenomena might also involve a very strong magnetic fields. Among them are gamma ray 
bursts and hypernovae It is supposed for the latter that a jet is somehow produced 
in these objects and the dipolar magnetic field is playing an essential role for that. On 
the other hand, some models for the gamma ray bursts are employing the magnetic fields 
to extract enormous energy of the phenomenon itself [Q. Furthermore, the evidence of 
generic asymmetry for collapse-driven supernovae has been accumulated [0, and the strong 
magnetic field might have some implications for the ordinary supernova if the magnetar 
as observed is the end product of the supernova explosion and the observed asymmetry is 
a dynamical consequence involving the strong magnetic field [Q. The fast proper motion 
of young pulsars might be explained by the combination of the strong magnetic fields and 
some processes such as neutrino oscillations, for example The explosion mechanism of 
the supernova will be changed substantially as well as nucleosynthesis therein. 

It is, therefore, not only of academic interest to consider the features of the strongly mag- 
netized vacuum. In particular, the quantum electrodynamical processes are most important, 
since those objects quoted above are mostly observed by electromagnetic waves, though the 
weak interactions are no less important [miO|- The study of the strongly magnetized vac- 
uum has a long history, though. Adler |Tl|, for example, gave detailed formulations for the 
polarization tensor and the photon splitting rate as well as some useful analytic expressions 
for limit cases (see also p|,|12|,|13[ and the references therein for many other contributions). 
Astrophysicists have used these approximate formulae for their model building |T^,|T3|. 

Those analytic expressions are approximate ones, though, valid for some limit cases such 
as the strong or weak magnetic field limits and the zero photon energy limit. It appears that 
we are lacking the complementary numerical evaluations of the polarization tensor for the 
intermediate values of magnetic field strength and/or photon energy. It is the purpose of the 
paper to fill this gap and give the interpolation formula based on the fitting to the numerical 
integrations. Our interest is, however, also directed to the dispersive relation of photon in the 
strongly magnetized plasma. In the last section we will extend the Schwinger's formulation 
and give the expression for the retarded polarization tensor for finite temperature plasmas. 



II. VACUUM POLARIZATION IN A STRONG MAGNETIC FIELD 



Using Schwinger's proper-time method |T^, we obtain the vacuum polarization tensor 
[see Fig. |l|) in a strong magnetic field expanded as. 
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(1) 



1=0 



where k = [uj, k±,0, k\\) is the energy-momentum 4- vector of photon, uj the photon energy, 
k± — a;sin^, k\\ — cucosO, theta the direction of 3-momentum with respect to the magnetic 
field, and B is the magnetic field. Gi is given by the following equations. 
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(6) 



where e denotes electron charge in MKS unit (e^ ~ 1/137), rUe is electron mass, and L{ 
B/Bc) is a dimensionless magnetic field normalized by the critical magnetic field [Be 
ml/e). f^" is expressed as 



/o — 9 



n 



k'^k'' 



(for/ = 1,2), 



with 



= k^Ff'^F^sk^ - k^'k^F^sF^'k,, 
K = k\ 



(7) 
(8) 



(9) 
(10) 

(11) 



where F*^^ denotes Maxwell stress 4-tensor and its dual tensor is defined by F**^^ = 

_1 /ii/7<5p 
2 c 'yd • 

It is easy for us to check that only Go contributes to the vacuum polarization tensor if 
the magnetic field is very week L <^ 1. In such a week field limit, we know the form of 
W{k) in usual QED as. 



W^{k,^)^G^\B=o!^\ 



(12) 
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where 



0\B=0 



(2vr) 



da dp 



2a 



exp 



-a + 



a(l-p^) 



(13) 



In Eq. (|1^) we should get rid of the divergence at a = and regularize it. Then we obtain 
the regularized form of the vacuum polarization tensor at L = 0, 



regU^^{k,0) = A{k)f^ 



with 



A{k) = regG'o|B=o 



(2vr) 



1 _ (1 - h{k)){4:ml + 2e) 
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where 



h{k) 



^-icot- y^i^ 



l_±5lcoth-M^/^ 



(for 0<P < Ami) 
(for A;2 < 0). 



(14) 



(15) 



(16) 



Thus, to obtain the regularized form of the polarization tensor in a strong magnetic field 
(L > 1), we have only to substitute Go with 



regGo = Go - Go|b=o + A{k). 



(17) 



III. REFRACTIVE INDICES IN A STRONG MAGNETIC FIELD 

As we mentioned in the previous sections, the refractive indices of photon would deviate 
from unity in a strong magnetic field because the vacuum polarization is influenced by the 
magnetic field and the dispersion relation is changed. The refractive indices is defined from 
the dispersion relation as, 



\k 



2 



/i' = (18) 



a;2 



where fc is a spatial 3-vector of k. To obtain the dispersion relation in a strong magnetic 
field, we consider the wave equation of photon, 

- k^k" + 47r regn'^'^(fc)] A^{k) = 0. (19) 

In this equation the prefactor of n^'^(/c), "47r", originates in MKS unit. The determinant of 
the matrix should be zero so that Eq. (|19D could have a nontrivial solution. When we choose 
the radiation gauge = (0, A), we get a quadratic equation and we obtain two solutions, 

1 + Xo-sm 6x1 
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^^-l + Xo + cos2^X2' ^ ^ 



where 



(22) 
(23) 

(24) 



Xo = 


47rregGo 








Xi = 


kl ' 




A7rG2 


X2 = 


kl ■ 



Here /i^ corresponds to the eigen vector (0,1,0) and to 
((1 + Xo + X2)cos^, 0, -(l + Xo)sin^). 

Since each Xi depends on jj? through k and Gj, Eqs. (^) and (|2T|) are imphcit equations 



for /i . In addition, in the hteratures, e.g. |T2[, they gave only the integral form and some 



limit cases. Hence we must solve the equations numerically to obtain values in the entire 
parameter space. In Fig. ^ we plot — xo as a function of L {=B /Be) in a low energy limit 
(a;^ = 10~^mg). From the plot, we find that the magnitude of Xo increases as L increases. 
Thus the contribution from xo to the refractive indices in Eqs. (^) would not be negligible 
in an extremely strong magnetic field. This behavior agrees with our analytical estimation 
that Xo oc — log(L) in a strong magnetic field limit {L oo) and contradicts with statement 
by Melrose and Stoneham |jl2[ that xo oc exp[— L] for this limit. This difference, however, is 



substantial only for extremely large magnetic fields. Nonetheless, we did not drop this term 
in estimating the correct refractive indices. In Fig. ^ we also plot Xi and X2 as a function of 
L. Xi approaches the limit value (~ 7.7 x 10^^) for L ^ 1. This feature is consistent with 
our analytical estimation and again disagrees with Melrose and Stoneham [jl2[ who claimed 
that Xi oc exp[— L] for this limit. As for X2, it is found that X2 is linearly proportional 
to L in the strong magnetic field and the weak energy limit. It agrees exactly with our 
analytical estimation that X2 — (e^/37r) L under the condition that L 1, a;^ — /cy <^ 4mg, 
and kl < 2mlL 

In Fig. ^ we plot the obtained refractive indices as a function of L. The solid and dashed 
lines represent ^\ and respectively. It is easy for us to understand the behavior of /i^s 



in a strong magnetic field. From Eq. ( P0|) and ( |2lD , we find that 

/i?~l + 7.7x lO-Sin^e, (25) 
/i^~l/cos^^, (26) 

for L ^ 1 in the case of the weak energy and as long as Xo is much smaller than unity. The 
photon-energy dependence of /x^ is shown in Fig. ^ In this plot we find that only near the 
threshold {uo'^ ~ Ami) ^^e energy dependence of /i^ becomes important, /if is insensitive to 
the photon energy. In Fig. |^ we plot function of cos ^ in a strong magnetic field 

(L = 10^). It is clear that /i^ is proportional to 1/cos^ 6. 
Here we give the fitting formula of 

/i2 = Ci(tanh[C2(x-C3)] + !) + !, (27) 

with 
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X = logio(S/Sc), 



(28) 



and 

r Ci = tan2^/2, 

I C2 = 1.15 - 7.07 X 10-3 (cuVm2)'-'° /V^^, (29) 

i C3 = 3.11 - logio (cos^ ^) - 1-84 X 10-2 (iuyml) / Vcos^ 9. 

This fitting formula reproduces tlie numerical results within the error of less than 10 % for 
a wide parameter range (0.1 < cos^ 6^ < 1, u^/ml < 4, and < B/ ^ 10^°). In the large 
L limit especially, Eq. ( pTj) approaches the value of the analytical estimation (—1/ cos^^). 
In addition, the low energy limit of the fitting formula (a;^ 0) agrees with the numerical 
estimations very well within less than 1 %. 



IV. RETARDED POLARIZATION TENSOR IN FINITE TEMPERATURE 

PLASMAS 

In the previous sections, we have considered the vacuum polarization tensor. It is also 
our concern to calculate the polarization tensor for plasmas with finite temperatures. We 
will extend the previous formulation to the finite density and temperature case. We will 
rely on the real time formalism of the finite density and temperature field theory and ob- 
tain the expression for the retarded polarization tensor. Recently, some authors |17] gave 



the expression of the chronological polarization tensor on a similar footing. Although two 
polarizations are related to each other, the retarded one allows more direct physical interpre- 
tation. Moreover, the retarded polarization tensor should be obtained by analytical extension 
of the imaginary time polarization tensor which has been given by the author ||18[ . 

It is known from the finite density and temperature field theory that the retarded polar- 
ization tensor is expressed as 

m^^ix, x') = Tr [YGc{x, x')YGr{x', x) + YGa{x, x')YGc{x', x)] . (30) 

In the above equation, the subscripts r, a, c denote the retarded, advanced and Keldysh 
components of Green function, respectively. For the stationary system, we can in general 
assume that the vector potential is time independent. In this case the polarization tensor 
Ilr{x, x') is a function of the time difference t — t' alone. It is possible that it depends on the 
spatial coordinates x and x' separately. Fourier transforming the polarization tensor with 
respect to t — t', we obtain 

illripo) = / ^ Tr [rG,{qo + Po)YGriqo) + 7^G,(go + Po)YGciqo)] ■ (31) 

In the above equation, the tilde means the Fourier component and the spatial coordinates 
are dropped for simplicity. Noting that the retarded and advanced Green functions for the 
finite density and temperature are identical to the counterparts for vacuum, we obtain 

1 f , G^''"(cj) - %^(o;) 



GrAPo) = 77- / , (32) 
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where the upper and lower signs correspond to the retarded and advanced Green functions, 
respectively. The subscripts F and F stand for the chronological and antichronological Green 
functions, respectively, for magnetized vacuum, which are calculated by the Schwinger's 
proper time method as shown above. It is also known that the Keldysh component of 
Green function can be obtained from the retarded and advanced Green functions and the 
distribution function by the following relation: 



Ge(po) = a(Po)[l - 2/(po)] - [1 - 2/(po)]G.(po) 



(33) 



Combining this with Eq. (|32|) , we can obtain the Keldysh component of Green function also 
from the chronological and anti-chronological Green functions for vacuum. 
Putting Eqs. (|32D, (|33D into Eq. (pTD and using the relations 

(34) 
(35) 



GMo) = Qiqo)GF{qo) - e(-go)G^(go) 
Gaiqo) = eiqo)Gp{qo) - e(-go)Gf (go) 



we finally obtain 
iflripo) = 



^[1 - 2/(go +Po)] X Tr [e(go)7'^G'^(go + Po)7'^G^(go) 

- Q{qo)l^Gp{qo+Po)YGF{qo) 

- ei-qo)YGF{qo+Po)YGpiqo) 
+ e{-qo)YG^{qo+Po)YG^{qo) 

^[1 - 2/(go -Po)] X Tr [e(-go)7'^G^(go)7'^G'^(go - Po) 

- e{-qo)rGF{qo)YG^{qo-po) 

- Qiqo)rGpiqo)YGFiqo-Po) 
+ Q{qo)YGp{qo)YGpiqo-po) 



(36) 



One easily recognizes that the structure of the integrand is quite similar to the vacuum 
polarization tensor apart from the integration over go and the various combinations of Gp 
and Gp. This enables us to simplify the integrand along the same line as for the vacuum 
case. 

Denoting as ^'^11^ the contribution from the terms containing the product GpGp and 
similarly for the other contributions, we calculate separately those terms. The antichrono- 
logical Green function is obtained by taking the integration region [0, — cxd] instead of 
[0, oo] in the Schwinger's proper time formalism : 



1 

i Jo 



dsexp -i (m'^ - (7^!^)^) s (7^!^ + m 



(37) 



Following Stoneham, we can simplify this equation. Assuming that the vector potential 
is time independent and Aq = 0, we can Fourier transform Gp and Gp with respect to 
t — t' . Plugging them into the definition of the polarization tensor, one sees that the gauge 
dependent terms cancel out just like the vacuum case and the polarization tensor becomes 
the function of the difference of the spatial coordinates, which then makes it possible for us 
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to take the Fourier transformation with respect to the spatial coordinates. We finally obtain 
for the contribution from the terms with the product GpGp 



dh 







X 



X exp 
X c/^^ 



oo 47rm 
^ da 
a 

cos l3 



FF{ko,Po\ 



dp exp 



.01 




cos a fc^ 



2L sin a 



exp 



exp 



.a 



—I- 



0^ k: 



a + P {ko +po)^ a- P kl 



2L 



+ 



2L 



Here FF{ko, po) is an abbreviation for the following function, 

FF{ko,Po) = e{ko) [1 - 2f{ko+po)] - e{-ko-po) [1 - 2f{ko)] , 
= e{ko) eiko + Po) [l-2Mko+Po)] 
+ e{-ko) ei-ko - Po) [1 - 2/e+(|A;o|)] 
- e(A;o) Qi-ko-po) [1 - 2/,+ (|A;o +Po|) + 1 - 2/e(A;o)] 



(38) 
(39) 

(40) 



where fe and /e+ are Fermi-Dirac distribution functions for electron and positron, respec- 
tively. Except for the integral over the distribution functions, the resemblance of the Eq. (pHf ) 
to the vacuum counter part is clear. The remaining factor d^" , which is symmetric with re- 
spective to the superscripts, is given as 



d 



'01 



1 cos P — COS a 

2 sin a 



+ 



j02 



1 cos P — COS a 

2 sin a 



+ 



cot 



cot 



cot 



cot 



'a — P'^ 
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'a — P^ 
'a + P^ 



cot a 



cot a 



cot a 



cot a 



ko+Po 
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h 
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ko+Po 



cot 



cot 



a 



^ a + P " 
P\ k^ 
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m 



cot 



'a + p' 



kx ^ ky 
m m 



m 



(41) 
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m 



m 



cot 



'a — P^ 



ky _|_ kx 
m m 



d'^ 
d^' 



1 k^ 
— cot a — 

2 m 



2A;o , Po 



+ 
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sin P 



sma 
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ko + Po ko 
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- P^ kl 



1 cos P 
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1 cos P 
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1 - 
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P tan/3 
a tan a 
P tan/3 
a tana 

ko + Po ko 
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m m 
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L 



4a2 



+ 



cos a — cos P kx k 
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sin^ a 



m m 



(42) 
(43) 
(44) 
(45) 
(46) 



m 



m 



1 + — ^ + 
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P^ kl 
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sin^ a 



Aa^ m 
1 cos P 
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sin^ a 
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(47) 
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Jll 



j22 
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sm a 



cos P 



sm a 



2a 



4a2 



cos /5 — cos a 



2 sin'^ a 



^0 + Po ko 
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2 sin^ a 
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(48) 



(49) 



j33 



cot a 



ko+Poko L a 

^^ + 7r^^ — 

mm la 4a^ 



15- kl 



m^ 



sin^ a 



1 cos /3 — cos a I kl k 



sin^ a 



m^ 



m^ 



(50) 



In the above equation, the integral region for /? can be changed from [—a, a] to [0, a] under 
the recognition that we take only the even part of the integrand with respect to (3. 

The other contributions to 11,. with different combinations of Gp and Gp are obtained 
in the same way. It turns out that the resultant equations are obtained from Eq. ( p8| ) with 
a change of the integral region and a replacement of the distribution functions as shown 
below: 



^^U'^'':FF{ko,Po)^FF{ko,Po) 
^^nr :FF(fco,j9o)^FF(fco,Po) 



FF 



integral region 
integral region 
integral region 



da 

-oo Ja 
oo 

da 

-oo J — CXD 



da I dp 



dp 
'dp 



(51) 
(52) 
(53) 



In the above equations, the phase should be taken as ^/—a = —i\/\a\, and the factors 
involving distribution functions are given as 



FF{ko,Po) = -Q{-ko) [1 - 2f{ko+Po)] - Q{-ko-po) [1 - 2f{ko)] 
FFiko,po) = -eiko) [l-2fiko + po)]-e{ko+po) [1 - 2/(A;o)] 
FFiko,Po) = ei-ko) [1 - 2fiko + Po)] + ei-ko-po) [1 - 2/(A;o)] 

Thus the polarization tensor is give by the sum of these terms: 



r 



(54) 
(55) 
(56) 



(57) 



Although the final form is similar to the vacuum counter part, the numerical evaluation 
is very difficult. This is so for the same reason as for the numerical evaluation of the 
vacuum polarization tensor above the threshold of the pair creation. This prevents us from 
performing a Wick rotation for the a and P integrations in Eq. (|38D . These should be the 
next step. 
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V. CONCLUSION 



We have numerically calculated the vacuum polarization tensor for various strengths 
of the background magnetic field and photon energies below the threshold of the pair cre- 
ation. We also varied the propagation directions of photon. We have obtained the fitting 
formula which reproduces the numerical results within 10% of error. Wc have also presented 
the expression of the retarded polarization tensor for the finite density and temperature, 
simplifying some integrations. 

Our final goal is to calculate not only the polarization tensor but also other physical 
quantities under strong magnetic fields such as the photon splitting rates and equation of 
states for a wide range of parameters. This paper is the first step for this project. 
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FIGURES 



wwv^yv\A/^ 




WWVVVW^ 



FIG. 1. Vacuum polarization in a strong magnetic field. The double solid line denotes the 
electron propagator which includes all contributions of the vertices from the magnetic field. 
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FIG. 2. Plots of — XOj Xi and X2 as a function of B/Bc- Here we adapt = 10~^me and 
cos^ = 0.5. We find that the magnitude of xo increases as B/Bc increases (oc \og{B / Be)), Xi 
reaches the hmit value (~ 7.7 xlO~^), and X2 — f^-^Z-^c for B/Bc > 1- 
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FIG. 3. Plot of the refractive indices as a function of B / Be- The soHd hne represents and 
the dashed hne represents ^j^- The left (right) dashed line is the case for u;^ = 4mg (10~^mg). Here 
we adapt = 10~^mg and cos^ = 0.5. We find that reaches the limit value (1/ cos^ = 2) in 
the strong B limit. 
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FIG. 4. Plot of the refractive indices as a function of lo^ /ml. The sohd line represents /if and 
the dashed line represents /i^- Here we adapt B/B^ = 10^ and cos^ 9 = 0.5. 
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FIG. 5. Plot of the refractive indices as a function of cos^O. The sohd hne represents /if and 
the dashed hne represents /x^- Here we adapt B/Bf. = 10^ and uj"^ = 10~^mg. We find that ^2 
scales 1/ cos^ in such a strong B limit. 
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